PASJ: Publ. Astron. Soc. Japan , 1-??, — 
© 2009. Astronomical Society of Japan. 



Suzaku Observation of the Metallicity 
in the Interstellar Medium of NGC 4258 
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Abstract 

The Suzaku X-ray satellite observed the nearby spiral galaxy NGC 4258 for a total good exposure time 
of 100 ks. We present an analysis of the Suzaku XIS data, in which we confirm that the 0.5-2 keV spectra 
of the interstellar medium (ISM) arc well-represented by a two-temperature model. The cool and hot ISM 
temperatures are 0.23tQ 02 and 0.59 ±0.01 keV, respectively. Suzaku's excellent spectral sensitivity enables 
us to measure the metal abundances of O, Ne, Mg, Si and Fe of the ISM for the first time. The resultant 
abundance pattern of O, Mg, Si, and Fe is consistent with that of the new solar abundance table of Loddcrs 
(2003), rather than Anders & Grevesse (1989). This suggests that the metal enrichment processes of NGC 
4258 and of our Galaxy are similar. 
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1. Introduction 

Metal abundances in the hot X-ray emitting interstellar 
medium (ISM) are important for understanding the star 
formation history and evolution of galaxies. A large frac- 
tion of metals in the ISM are synthesized by type la and 
type II supernovae (SNe la/II). O and Mg are predomi- 
nantly synthesized by SNe II, while Fe is mainly produced 
by SNe la. Therefore, the abundance ratios provide use- 
ful information on the contribution of both types of SNe 
enriching the metals. 

Compared to starburst galaxies, normal spiral galax- 
ies have lower star formation rates (SFRs), and have 
lower X-ray surface brightness, because the integrated soft 
X-ray luminosities of spiral galaxies correlate with SFR 
(Tiillmann et al. 2006). In addition, these galaxies some- 
times exhibit strong X-ray emission from the central active 
galactic nuclei. As a result, it has remained rather diffi- 
cult to constrain the abundances of their hot ISM over a 
wide range of species from O to Fe. 

In order to determine the metal abundances of O, Ne, S, 
and Ar from the emission lines in H II regions, several opti- 
cal observations have been performed (e.g., Perez-Montero 
et al. 2007, Crockett et al. 2006). In particular, O emis- 
sion lines, a good tracer of SNe II, are the brightest in the 
optical band. However, determination of the metal abun- 
dances, including O, has large uncertainty due to dust 
grains, the composition of which is difficult to ascertain. 
In addition, no significant Fe lines, a good tracer of SN 
la, exist in the optical band. 



The Einstein satellite first observed the thermal X- 
ray emission from spiral galaxies (Fabbiano 1989), and 
ROSAT and ASCA subsequently measured the tem- 
perature and Fe ISM abundances (Roberts & Warwick 
2001; Petre et al. 1994). Chandra and XMM-Newton have 
enabled us to study the spatial distribution of metals by 
utilizing their large effective area and high angular reso- 
lution (Schlegel et al. 2003; Yang et al. 2007). However, 
the O and Mg abundance measurements are difficult due 
to the intrinsic instrumental background and asymmetric 
energy response in the low energy band below ~ 1 keV. 
Because Suzaku XIS (Koyama et al. 2007) has higher spec- 
tral sensitivity below ~ 1 keV, and a lower and more stable 
background level compared to Chandra ACTS and XMM- 
Newton EPIC, the accuracy of the determination of O 
and Mg abundances is improved. Recent Suzaku obser- 
vations have revealed the abundance profiles of the ISM 
in X-ray bright elliptical and starburst galaxies in detail 
(Matsushita et al. 2007; Tsuru et al. 2007; Tawara et al. 
2008; Yamasaki et al. 2009; Hayashi 2009). 

NGC 4258 (M 106) is a nearby SABbc spiral galaxy 
with a high inclination angle of 72° (TuUy 1988). Its dis- 
tance is accurately measured to be 7.2 Mpc (Herrnstein et 
al. 1999), where 1' corresponds to 2.1 kpc. Makishima et 
al. (1994) first discovered the obscured low-luminosity ac- 
tive galactic nucleus (LLAGN) in NGC 4258 with ASCA, 
and also showed that its 1-10 keV spectrum is composed of 
contributions from thermal plasmas, integrated low mass 
X-ray binaries (LMXBs), and the nucleus. A ROSAT ob- 
servation detected two thermal components with kT ~ 0.2 
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Fig. 1. The X-ray contour map (orange) in linear scale from 0.5-2 keV (left) and 2-10 keV (right), overlaid on an optical image 
taken by DSS. The contour scales range linear by (0.17 - 5)xl0-^ and 0.3 - 1.3 counts sec ^ arcmin ^, respectively. The observed 
XISO, 1, 2, and 3 images were added on the sky coordinate, and smoothed with a = 12 pixel Gaussian profile. Estimated contributions 
of the cosmic X-ray background (CXB) and the instrumental non X-ray background (NXB) were subtracted, although vignetting 
was not corrected. The white square corresponds to the field of view of the Suzaku XIS. The region where energy spectra were 
extracted is indicated by a white dashed ellipse. Three point sources (white dashed circles) were excluded from the background 
energy spectra. 



and - 0.5 kcV in the ISM (Voglcr & Pictsch 1999), which 
were confirmed with Chandra and XMM-Newton observa- 
tions (Yang et al. 2007). The 2-10 keV flux from the nu- 
cleus was reported to vary on time scales of hours to years 
(Reynolds et al. 2000; Tcrashima et al. 2002; Fruscione et 
al. 2005). A recent Suzaku observation has revealed that 
the circum-nuclear matter forms a slim obscuring struc- 
ture with solid angle much smaller than typical thick tori 
(Yamada et al. 2009; Reynolds et al. 2008). The galaxy 
has a curious jet-like feature, what is called an "anomalous 
arms", emitting soft X-rays (e.g., Cecil et al. 1992, Wilkes 
et al. 1995; Cecil et al. 1995). In addition, the hot spots 
associated with the galactic radio jet have been detected 
at a distance of around 1' from the nucleus in X-rays by 
Chandra (Yang et al. 2007). Aside from these anoma- 
lous arms, which seem to be related to the nuclear activ- 
ity, NGC 4258 can be regarded as being similar to our 
Galaxy in terms of the low activity of the central black 
hole. The fact that the emission from the nucleus is deeply 
absorbed with a thick circum-nuclear matter (Makishima 
et al. 1994) is helpful to study the thermal emission. In 
this paper, we report the detailed analysis of NGC 4258 
with Suzaku, focusing on the thermal emission of the ISM 
to constrain the metal abundances. 

Unless noted otherwise, the solar abundance table used 
is taken from given by Anders & Grevesse (1989), and the 
quoted errors are for a 90% confidence interval for a single 
interesting parameter. 



2. Observation and Data Reduction 

Suzaku carried out an observation of NGC 4258 in 2006 
June. We analyzed only the XIS data in this paper, al- 
though Suzaku observed the object with both XIS and 
HXD (Takahashi ct al. 2007), centered on the HXD field 
of view, and the HXD results from this observation on 
the nucleus have already been published by Yamada et al. 
(2009). The XIS consists of four sets of X-ray CCD camera 
system (XISO, 1, 2, and 3). XISl has a back-illuminated 
(BI) sensor, while XISO, 2, and 3 have front-illuminated 
(FI) sensors. The averaged pointing direction of the XIS 
is centered on (R.A., Dcc.)=(12'^18'"50?8, +47°14'15"4). 
The XIS was operated in the normal clocking mode (8 
s exposure per frame), with the standard 5x5 and 3 
X 3 editing mode. We processed the XIS data by the 
"xispi" and "makepi" ftools task and the CALDB files 
of 2008-01-31 version. Then, the XIS data were cleaned 
by assuming thresholds on the Earth elevation angle of 
> 5° and the Day-earth elevation angle of > 20° . We also 
discarded data during the time since the spacecraft exit 
from the south Atlantic anomaly was unless than 436 sec. 
We created a light curve of each sensor over 0.5-2 keV 
binned 540 sec to exclude periods of an anomalous event 
rate greater or less than ±3ct around the mean. After the 
above screening, the remaining good exposures were 99.9 
and 99.8 ksec for FIs and BI, respectively. Event screening 
with cut-off rigidity was not performed. Wc also excluded 
the RAWY<119 region on the XIS2 and 3 due to an er- 
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roneous dark frame during this observation. 

The spectral analysis was performed with HEAsoft ver- 
sion 6.5 and XSPEC 12.4. In order to subtract the non- 
X-ray background (NXB), wc employed the dark-Earth 
database by the "xisnxbgen" ftools task (Tawa et al. 
2008). We generated two different Ancillary Response 
Files (ARFs) for the spectrum of each region, which as- 
sumed uniform sky emission and the observed XISl image 
by the "xissimarfgen" ftools task (Ishisaki et al. 2007). We 
also included the effect of contaminations on the optical 
blocking filters of the XIS in the ARFs. 

3. Analysis and Results 

3.1. Strategy of Spectral Fitting 

Figure 1 shows the observed X-ray contours in the 0.5-2 
keV and 2-10 keV ranges, overlaid on the optical image 
from the Digitized Sky Survey (DSS). We extracted spec- 
tra from the following two regions. 

1. NGC 4258 components: an ellipse region, with semi- 
major and semiminor axes of 8.3' and 3.6'. 

2. Background: the entire XIS field of view was used, 
after excluding the above NGC 4258 component and 
three additional 1.5' radius circular point source re- 
gions mentioned below. 

Three X-ray point-like sources were detected in 0.5-2 
keV as shown in figure 1. Sources #1, 2, and 3 were 
identified as a quasar Q 1218-1-472 (Burbidge 1995), an 
Ultra Luminous X-ray source (Arp et al. 2004), and an 
irregular galaxy UGC 7356 (Thuan & Seitzer 1979), re- 
spectively. Although these sources were not so bright as 
to affect the spectra of NGC 4258, these were not negli- 
gible with respect to the estimation of the extra-galactic 
cosmic X-ray background (CXB) and the Galactic emis- 
sions. In addition, three calibration sources, which have 
emission peaks at 5.9 keV, are located at the corners of 
the XIS. We included these regions for improvement in 
photon statistics because our target spectra are extracted 
only over the energy range below 5 keV. 

The 0.5-2 keV contours of NGC 4258, which mainly 
represent the ISM emission, are clearly extended as com- 
pared to those in 2-10 keV, which are dominated by the 
LLAGN emission. Therefore, the soft X-ray signal is con- 
sidered to arise mainly from an extended hot halo. In fact, 
Yang et al. (2007) showed that the 0.2-10 keV spectra ob- 
served with XMM-Newton were represented by the sum 
of two thermal components of 0.60 and 0.22 keV, plus the 
emission from LMXBs and the obscured LLAGN. The re- 
sultant parameters of these components imply that, even 
in the nuclear region within 3', the diffuse thermal ISM 
emission dominate the total flux below ^ 2 keV. Thus, we 
also analyzed the XIS spectra below 2 keV to investigate 
the ISM properties. 

3.2. Estimation of the background spectra 

In order to estimate the CXB and Galactic emission 
background components we first fitted the spectra of the 
background region. We assumed a power-law model for 



the CXB component, and a two temperature model for the 
Galaxy to represent emission from local hot bubble (LHB) 
and Milky Way halo (MWH). We tested the following two 
models: (i) phabs x power-law -I- apecMWH + apccLHB, 
and (ii) phabs x (power-law -I- apecMwn) + apecLHB, 
where the apec models have a metal abundance fixed at 
the solar value and zero redshift. The absorption column 
density was also fixed to the Galactic value in the direc- 
tion of NGC 4258, which is 1.2x102" cm'^. The difference 
between the two models is whether the MWH component 
is absorbed by the Galactic column or not (Yamasaki et 
al. 2009). We also fixed the temperature of the LHB and 
MWH emission to be 0.1 and 0.3 keV, respectively, af- 
ter previous studies (Lumb et al. 2002; Yamasaki et al. 
2009; Komiyama et al. 2009). Due to poor photon statis- 
tics, leaving either temperature free did not provide use- 
ful constraints. The spectra from the BI and FI CCDs 
were fitted simultaneously in the 0.5-5.0 keV range. The 
results of these fits are summarized in table 1. The nor- 
malization of the power-law model was found to be higher 
(~19%) than the previously reported value (e.g., Kushino 
et al. 2002). This may be accounted for as due to the 
CXB brightness fluctuation, or to possible contamination 
by point sources in NGC 4258. Although the results do 
not differ significantly between the two models, model (ii) 
gave a marg inally better fit statistic, x^/d.o.f. = 558/426. 

In order to estimate the Galactic emission, we fitted 
the OVII and OVIII lines with Gaussian and a power- 
law in the 0.5-0.7 keV range. The resultant intensities 
of the O VII and O VIII without correcting for absorption 
are 3.8±1.1 and 1.1±0.9 photons cm"^ s~^ sr~^, respec- 
tively. These are consistent with previously reported val- 
ues (McCammon et al. 2002; Sato et al. 2007; Yamasaki 
et al. 2009). Hereafter, model (ii) is utilized as the back- 
ground for the spectral fitting, unless otherwise stated. 

3.3. The galaxy spectra of NGC 4258 

We fitted the galaxy spectra with the model: 
phabsG X (vapeciT,2T,or3T + bremss) -I- phabs agnX 
power-law AGN + phabsGX (power-lawcxB + apecMwn) 
+ apccLHB , where the last term represents the background 
model to be explained later. In the model, phabsG means 
the Galactic absorption in the direction of NGC 4258, 
fixed at Nu = 1.2 x 10^° cm"^. The term (phabsAGNX 
power-law agn) shows the LLAGN contribution; its ab- 
sorbing column was fixed at 1.07x10^° cm^^, and the 
power-law AGN slope and normalization at F = 1.86 and 
Norm = 4.22 x 10"'"^ photons keV~^ cm"^ s~^ at 1 keV, 
both after Yamada et al. (2009) . This formalism assumes 
that the absorption within NGC 4258 follows the same 
abundance ratios as the Galactic ones employed in the 
phabs model. This assumption is justified a posteriori by 
our result to be obtained below, that NGC 4258 and our 
Galaxy have rather similar abundance patterns. Although 
the intensity of the LLAGN component is much lower 
than that of the ISM components below 2 keV (subsection 
3.1), we take this term into account to decide accurately 
the emission of the Si-Ka line around 1.8 keV where the 
LLAGN emission may affect the spectral fits. The term 
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Table 1. The best-fit parameters of the apec component + pawer-law model.* 



Parameters 




(i) 


(ii) 


TcxB 

NormcxB ^ 

/Vw 
J V fi 




1.38±0.05 
1 16+005 
^■^"J-o.oe 
1 2x1 n^o (fi^] 


1.37±0.05 

1 1 c+0.05 
^■^^-0.06 

1 2x1 0^0 (fix) 


Abundance 
NorniMWH * 


(keV) 
(solar) 


0.3 (fix) 
1 (fix) 
0.64±0.06 


0.3 (fix) 
1 (fix) 
0.68±0.06 


Abundance 
A^ormLHB ^ 


(keV) 
(solar) 


0.1 (fix) 
1 (fix) 

1 74+0.36 


0.1 (fix) 
1 (fix) 
1.80t°;f, 


xVd.o.i. 




562/426 


558/426 



* The apec components for the spectra in the background region of NGC 4258 with absorbed or non-absorbed MWH component for the 
Galactic emission, and a power-law model for CXB. 

t Normalization of the power-law component divided by the solid angle, f2", assumed in the uniform-sky ARF calculation (20' radius), in 
units of 10""^ r2" photons keV~^ cm~^ arcmin"^ at 1 keV. 

t Normalization of the apec components divided by the solid angle same as the normalization of apec, Norm = J ncTijidV / (47r(l + 
2)^D^) /fl" xlO~^^ cm~* arcmin"^, where Da is the angular distance to the source. 



(phabscx (power-lawcxB + apecMwn) + apecLHB) rep- 
resents the background component as examined in sub- 
section 3.2. 

The ISM emission of NGC 4258 was modeled with one, 
two, or three-temperature models, as indicated by the sub- 
scripts IT, 2T, or 3T, respectively, employing the vapec 
code (Smith et al. 2001). The abundances of He, C, and 
N were fixed to the solar value. We also divided the 
other metals into five groups as O, Ne, (Mg & Al), (Si, 
S, Ar, Ca), and (Fe & Ni), based on the metal synthesis 
mechanism of SNe, and allowed them to vary. The abun- 
dances were constrained to be common all temperature 
components. The bremss model, with kT= 10 keV, repre- 
sents the integrated LMXB component (Makishima et al. 
1994; Yamada et al. 2009). In order to constrain the back- 
ground component contained in the above fitting model, 
we simultaneously fitted the source and background spec- 
tra, over the 0.5-2 keV and 0.5-5 kcV regions, respec- 
tively. When fitting the background spectra, the normal- 
izations of vapec, bremss, and power-law agn were all fixed 
to be 0. 

The results of these fits are shown in figure 2, and the 
derived parameters are summarized in table 2. Thus, 
several emission lines arc seen around 0.5-0.6 keV, 0.6- 
0.7 keV, ~ 1.3 keV, and ~ 1.8 keV, and are identified 
with those of O VII, O VIII, MgXI and Si XIII, respectively. 
Furthermore, emission lines around 0.7-1 keV correspond 
to the Fe-L complex, as well as to K-lines from Ne IX and 
NeX. The fits are not formally acceptable because of the 
big residuals around 2 keV owing to calibration issues of 
Si edges, which are known to be present in all the XIS 
CCDs (Koyama et al. 2007). However, these results are 
still useful to assess whether the temperatures and metal 
abundances in the ISM were reasonably determined or 
not. The one-temperature model results in kT = 0.38 keV. 
The fit statistics shown in table 2 clearly favor the two- 



temperature model, which gives the two temperatures as 
0.23 and 0.59 keV, with corresponding 0.2-2.4 keV un- 
absorbed fluxes of 1.0x10"^^ ^j^^ 2.5x10"^^ ^g cm"^ 
s^\ respectively. The 0.1-2.4 keV total unabsorbcd ISM 
luminosity is 2.2x10^" erg s^^, which is consistent with 
that of the previous XMM-Newton results (Yang et al. 
2007). The three-temperature model improved the fit 
statistics only slightly, and yielded the temperatures as 
0.13, 0.36, and 0.62 keV. The implied abundance ratios 
did not change, as shown in table 2. 

In order to examine abundance ratios rather than ab- 
solute values, we calculated confidence contours between 
the abundance of metals (O, Ne, Mg and Si) and that of 
Fe, using the two-temperature model for the ISM. The 
results are shown in figure 3, where we also indicate 90%- 
confidcncc abundance of these metals relative to Fe. 

3.4- Systematic Uncertainties in the Abundance Ratios 

In order to investigate how these results are sensitive 
to the assumed Galactic model, we let the temperatures 
of the MWH and LHB vary freely. Although the temper- 
atures of the Galactic components became 0.17 and 0.61 
keV, the ISM abundance ratios were little affected. When 
we change the normalizations of Galactic apecMWH and 
apecLHB by ±10% and ±20%, respectively, the metals to 
Fe abundance ratios of the ISM remained the same within 
^5%. The systematic uncertainty due to the assumption 
of the OBF contaminant is less than the statistical error 
in the NGC 4258 region. Thus, the uncertainties in the 
background and the OBF contaminant have little effect 
on the abundance ratios of the NGC 4258 ISM. 

Next, we investigate uncertainties in the absorption pa- 
rameter. Although we have carried out the fits with the 
absorption fixed to the Galactic value, Yang et al. (2007) 
have shown that absorption varies by a factor of a few, 
depending on location. When we fit the spectra with the 
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Fig. 2. The NXB-subtractcd XIS-FI (black) and XIS-BI (red) spectra of NGC 4258 (left column), and those of the background 
region, shown without removing the instrumental responses. Each spectrum from the top to bottom panel was fitted with a one, 
two, or three-temperature model for the ISM, respectively. The black and red lines show the best-fit model for the FI and BI, 
respectively. For simplicity, only the model components for BI spectra are shown. The magenta, orange, and purple lines are the 
ISM components, blue and cyan are the Galactic background emission by apecMWH and apecLHBi and green, gray, and yellow are 
the CXB, LMXB, and LLAGN components, respectively. The background components are common between the on-source and 
background spectra. 
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Table 2. Summary of the best-fit parameters of the NGC 4258 and background region. 



Parameters 



IT for ISM 2T for ISM 



TcXB 

NormcxB * 
kTuwH 
NoruiMWH ^ 

Normi^HB ^ 



kTn 

Normi'Y t 

kT2T 

Norm2T ^ 

NormyY ^ 

O 

Ne 

Mg, Al 
Si,S,Ar,Ca 
Fe, Ni 

LMXB flux* 
xVd.o.f. 



(keV) 
(keV) 

(keV) 

(keV) 

(keV) 

(solar) 
(solar) 
(solar) 
(solar) 
(solar) 



2T for ISM 
(LHB and MWH free) 



3T for ISM 



1.42 
1.21 
0.3 (fix) 
0.68 (fix) 
0.1 (fix) 
1.80 (fix) 

0.38±0.01 
2.18±0.04 



'?fi+"-i" 

1 20+°-^^ 

^■^"-0.20 

' ^-0.15 
1+0.54 
-0.39 

4q+" i2 

"■^^-0.08 



1.59: 



1.41±0.04 
1.21±0.04 

0.3 (fix) 
0.68 (fix) 

0.1 (fix) 
1.80 (fix) 

0.23±0.01 
0.681^;}^ 
0.56±0.01 

1 1 c:+0.06 
J^-J^3-0.42 



n 56+''-i^ 

^•'JD-0.08 
1 90+0. 62 
^■^°-0.20 
-1+0.47 
-0.15 



0.99^ 

1 1 1+0.50 
J^-^^-0.18 

n 64+''-2'^ 

^•"J^-o.og 



(ergcni-2 s'^) 9.2 xlO'^^ 7 9 ^IQ- 



13 



1333/740 993/738 



1.31] 
1.08^ 



I+U.U4 
-0.06 
j+0.03 
'-0.02 

^•"-■^-0.03 

n 09+0.03 

^■'J^-0.04 

0.17±0.01 

0.88i°:°^ 



0.24±0.01 

^•'^-0.05 

0.58±0.01 

-1 97+0.26 
^■'^ ' -0.03 





'J-O'J-0.04 
1 -19+0.06 
^^•-■^^-0.07 
+0.13 
0.11 
+0.13 
0.11 
+0.01 
-0.02 



0.87 
0.99 
0.57 



7.7 xlO 



-13 



946/734 



1.41±0.04 
1.21±0.04 
0.3 (fix) 
0.68 (fix) 
0.1 (fix) 
1.80 (fix) 

'^•-'^'3-0.06 

41+011 

^^•^-^-0.02 

36+0 06 

'^•«J"-0.05 

fS2+0 03 
'-'•0^-0.02 

62+0 06 

^•"^-0.02 

3+0.04 
-0.02 
1+0.24 
-0.32 

-L-^o_0.46 
-, 9r:+0.42 

0.84tojo 



0.68] 
0.91^ 



7.8 xlO 



-13 



982/736 



* Normahzation of the power-law component divided by the solid angle, Q' , assumed in the uniform-sky ARF calculation (20' radius), in 
units of 10""^ r2" photons keV~l cm~^ arcmin"^ 

t Normalization of the vapec component scaled with a factor of SOURCE_r,atio_r,eg / AREA, which is Norm 

-17 5 arcmin"-^, where Da is the angular distance to the source. 



SOURCE_RATIO_REG 



AREA 



cm " s ^ arcmm ^ at 1 keV. 
vapec component scaled with 
J UeniidV / [Air (1 + z)^ Dl] xlO"" cm 



Flux within the accumulated region between 0.5 and 2 keV. 



absorption allowed to vary freely, the resultant absorp- 
tion is constrained with an upper limit of A'h of 4.7 x lO^o 
cm^^. Thus, wc confirm that the effect of the absorp- 
tion is almost negligible because the resultant abundance 
change insignificantly. 

In order to examine the dependence on the LMXB com- 
ponent, we used power-law model instead of the bremss 
model in Makishima et al. (1994) and Yamada et al. 
(2009). However, the abundance ratios to Fe did not 
change significantly either. 

The derived Ne to Fe abundance ratio is higher by a 
factor of ~ 2 higher than the solar abundance. This is 
presumably because the Ne abundance is not reliably de- 
termined due to an overlap with the strong and complex 
Fe-L line emissions. 

Yang et al. (2007) have detected the hot spots of the 
jet located around 1' from the nucleus in X-rays, and 
also found a non-thermal component besides that from 
the AGN within a 20" region centered on the nucleus. 
In order to examine the contribution of the jet emis- 
sion, we tried to investigate the spatial distribution of the 
metal abundances. Because the above physical scales are 
much smaller than the nominal 3' point spread function of 
Suzaku (Serlemitsos et al. 2007), we are only able to place 
crude constraints. We first extracted spectra from two 



separate regions for comparison: (i) a circular aperture of 
radius 3' centered on the nucleus, and (ii) the remainder 
of the elliptical data-accumulation region from figure 1, 
excluding the central 3'. We fitted the two spectra with 
the two-tempcrature model. Unfortunately, poor photon 
statistics prevented the detection of any significant abun- 
dance differences between the two regions. Thus wc car- 
ried out a different test in which we added the following 
jet model from Yang et al. (2007) to the two-temperature 
model : phabsjot x power- lawjet. Here, phabsjet, the in- 
trinsic absorption, was fixed at 3.4x10^^ cm~^, and the 
photon index of power-lawjct was fixed at 1.74 (Yang et al. 
2007). Then, the 0.5-2 keV jet emission flux was rather 
tightly constrained to be < 4.5xl0~i'' erg cm^^ s^^, at 
a confidence limit of 90%. Because the other parameters, 
including the metal abundances, did not change signifi- 
cantly, wc conclude that our results are not affected by 
the jet emission. 

4. Discussion 

The present Suzaku observation has clearly revealed 
emission line features from the ISM in the spiral galaxy 
NGC 4258. We successfully measured the metal abun- 
dances of O, Ne, Mg, Si and Fe for the first time. Previous 
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Fig. 3. Plot of the confidence contours between the metal (O, Mg, S, Si) and Fe abundances, indicated by the two-tempcrature 
model for the ISM. The three contours represent 68%, 90%, and 99% confidence range from inner to outer, respectively. 



observations with XMM-Newton and Chandra with better 
angular resolutions successfully measured spatial temper- 
ature variations and the geometric structure of "anoma- 
lous arms". But, the metal abundances were not con- 
strained due to insufficient photon statistics (Yang et al. 
2007). Figure 4 shows our mctal-to-Fc ratios, which were 
derived from the two-parameter confidence contours in fig- 
ure 3. The ISM emission is here modeled with two tem- 
peratures. 

In figure 4, the abundance pattern of NGC 4258 is com- 
pared with those indicated by the solar abundance table of 
Anders & Grevesse (1989) and the new solar abundance 
table of Lodders (2003). The O abundance in Lodders 
(2003) was derived from solar photosphcric lines, consider- 
ing three-dimcntional hydrostatic model atmospheres and 
non-local thermodynamic equilibrium (Asplund 2005 and 
references therein). Thus, the abundance pattern of the 
ISM of NGC 4258 measured with Suzaku agrees better 
with that of Lodders (2003), rather than that of Anders 
& Grevesse (1989). 

The calculated SN II and SN la yields are also plotted 



in figure 4. Here, the SN II yields by Nomoto et al. (2006) 
refer to an average over the Salpeter initial mass function 
of stellar masses from 10 to 50 Mq, with a progenitor 
metallicity oi Z = 0.02. The SN la yields were taken from 
the W7 model Iwamoto et al. (1999). The abundance 
pattern of the ISM of NGC 4258 is located between those 
of SN II and SN la, and consistent with their mixture is. 

In order to investigate differences between spiral and 
starburst galaxies, we also plot in figure 4 the abundance 
pattern of the hot ISM in disk and halo regions of the 
starburst galaxy NGC 4631 (Yamasaki et al. 2009), and 
that of the "cap" region of the extreme starburst galaxy 
M 82 (Tsuru et al. 2007). The results of NGC 4631 disk 
are consistent with those of Lodders (2003) as is the case 
of NGC 4258, while the respective patterns of M 82 "cap" 
and the halo region of NGC 4631 are close to those ex- 
pected from SN II yields. 

The SFR of NGC 4258 derived from its far infrared lu- 
minosity is low, '^O.Ol Mq yr~^ (Wu & Cao 2006), com- 
pared with those of NGC 4631 and M 82, which arc 2.99 
and 9.39 Mq yr~^, respectively (Tiillmann et al. 2006). 
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Fig. 4. Number ratios of O, Ne, Mg, Si and S to Fe for the two and three-temperature model of the ISM in NGC 4258. Solid 
and dashed lines indicate the number ratios of metals to Fe due to Lodders (2003) and Anders & Grevesse (1989), respectively. 
Abundance patterns of NGC 4631 disk, halo (Yamasaki et al. 2009), and M 82 "cap" (Tsuru et al. 2007) also are shown. Dot-Dashed 
and dotted lines represent the number ratios of metals to Fe for the SN II and SN la products (Iwamoto et al. 1999; Nomoto et al. 
2006). 



The 0.5-2 keV luminosity of the LMXB component of 
NGC 4258, 3.5 x 10^^ erg s~^, is consistent with the mass 
vs. X-ray luminosity relation for spiral galaxies (Gilfanov 
2004) . Thus, we can regard NGC 4258 as a "normal" spi- 
ral galaxy. Therefore, we may conclude that solar abun- 
dance pattern are common in normal spiral galaxies, in- 
cluding NGC 4258 and our Galaxy. The fact that the 
disk region of NGC 4631 also has a similar abundance 
pattern to normal spirals, suggests that the metallicity 
of the ISM after its starburst era may look quite similar. 
In contrast, in starburst galaxies, SN II products such 
as O effectively escape into the halo region as a result of 
the energetic explosions. Thus, observations of abundance 
patterns such as ours play a key role in investigating the 
processes of galaxy evolution and enrichment of the inter- 
galactic medium. 
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